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A new quasi-neutral model describing the Weibel instability of a high- current relativistic beam 
propagating through a resistive plasma is developed. It treats beam electrons as kinetic particles, and 
ambient plasma as a non-relativistic fluid. For a finite-temperature beam, a new class of negative 
energy magneto- sound waves is identified, which can possess negative energy. Their growth due 
to collisional dissipation in the cold return current destabilizes the beam-plasma system even for 
high beam temperatures. We perform detailed two- and three-dimensional particle-in-cell (PIC) 
simulations of the thermal beam and collisional plasma. It is shown that in three dimensions, 
the Weibel instability persists even for physically collisionless background plasma. The anomalous 
plasma resistivity is then caused by the two-stream instability. 
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The fast ignition fusion (FI) is a promising route to- 
wards the laser driven fusion. In the classic FI scheme 
a laser-generated relativistic electron beam with a 
few MeV per electron energy must propagate through 
overdense plasma to heat a hot spot in the core of a pre- 
compressed fusion fuel target. The current carried by 
these MeV electrons inside the plasma is much higher 
than the Alfven current limit / = (mc^/e)7 = 177/cA, 
where m is the electron mass, e is the electronic charge, 
and 7 is the Lorentz factor of the beam. Transporta- 
tion of this electron beam is not possible unless it is 
compensated by a return plasma current. However, this 
configuration is unstable and the current beam is sub- 
ject to the Weibel and the two-stream instabilities. The 
Weibel instability [2] is particularly responsible for the 
generation of very strong static magnetic fields (~ 100 
MG) i. It is one of the leading instabilities under rela- 
tivistic conditions and has been studied for a long time 
0,i,[i,[7, 8, 9, id El, El. Honrubia et al [11] have per- 
formed three-dimensional simulations of resistive beam 
filamentation corresponding to the full scale FI configu- 
ration. Three-dimensional magnetic structures generated 
due to the Weibel instability in a collisionless plasma 
have also been reported [12]. Recently, the evidence of 
Weibel-like dynamics and the resultant filamentation of 
electron beams have been reported experimentally [6] . It 
was proposed in Ref. [7] that this instability could be 
suppressed by the transverse beam temperature alone in 
a collisionless plasma. However, the instability persists 
in the presence of collisions in return plasma current no 
matter how high transverse beam temperature is. This 
regime of instability was termed as resistive beam insta- 
bility fi^. 

In this Letter, we develop a theoretical approach to 
the collisional Weibel instability in the framework of the 
quasi-neutrality assumption. For a finite-temperature 
beam, a new class of negative energy magneto-sound 



waves is identified. We derive conservation laws for the 
beam-plasma system and show that the energy of this 
system is not positively definite. Rather, it contains a 
negative term that allows for negative energy waves. Col- 
lisions in the background plasma current excite unstable 
magneto-sound waves in the system, which carry negative 
energy densities. Due to these waves the Weibel insta- 
bility is not suppressed even when the transverse beam 
temperature would be high enough to stabilize collision- 
less plasma. 

We present results of detailed 2D and 3D particle-in- 
cell (PIC) simulations on the relativistic electron beam 
transport in plasmas. The 2D geometry corresponds to 
a plane transverse to the beam propagation direction. In 
this geometry the Weibel instability is decoupled from 
the two-stream instability and we can study the effects 
of temperature and collisions on the Weibel instability 
systematically. The simulations results show that the 
Weibel instability cannot be suppressed by thermal ef- 
fects alone if collisions are present in the system. We 
also make 3D PIC simulations of the Weibel instability. 
The simulation results show that in the full 3D geome- 
try, the Weibel instability cannot be suppressed even in 
plasma free from binary collisions. We conjecture that 
the effective collisions leading to an anomalous resistiv- 
ity in the return current are provided by the turbulence 
emerging from the electrostatic beam instability. 

In the 2D simulations, we assume a very long electron 
beam propagating in the Z— direction and there is no 
dependence on the coordinate z. The beam and plasma 
densities are and Up respectively and the beam-plasma 
system is quasi- neutral i.e. nij -\- Up = no, where no is 
the background ion density. Initially, the beam current 
is completely neutralized by the return plasma current. 
However, as the Weibel instability develops the current 
neutrality is destroyed due to the filamentation. The 
strongest magnetic field is generated in the transverse 
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plane {x — y plane). This magnetic field generates an 
axial component of the electric field E^. The transverse 
components of the electric field can be easily obtained 



from the force equilibration E 



X B±/c = 0, where 



Bz is the return current velocity. To summa- 
rize, these are the dominant electric and magnetic fields 
of the beam plasma system: 



' c dt 



Bj_ = -e^xV±A^, = E± = -{vp^/c)V±A^, 

(1) 

where Az is the z component of the vector potential. The 
Bz component is small, but can be approximated using 
dtBz = — C^_LVpz X \/_\_Az) • Cz. The Az component is 
determined from the Ampere law 



Az 



47r 



{Jbz H~ Jpz) 1 



(2) 



where Jhz and Jpz are the current densities of the beam 
and plasma respectively. We discard the displacement 
current to ensure the quasi-neutrality. The axial equa- 
tion of motion for the plasma fiow and the transverse 
equations of motion for beam electrons are 
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where u is the collisional frequency of the ambient 
plasma, m and c are the electron mass and velocity of 
light in vacuum respectively, and the subscript j rep- 
resents the j^^ beam electron. For collisionless plasma, 
Eq. (|4]) is written as 
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where Lz is the system length in the z direction. The first 
term in the above expression represents the total beam 
electron energy. The third and fourth terms correspond 
to the total magnetic energy and the plasma kinetic en- 
ergy. The fourth term slightly overestimates the energy 
because the actual plasma density = no — rib is slightly 
smaller. The second term corrects this overcount: the ex- 
cess energy is subtracted from the electron beam energy. 
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The second term in the RHS of Eq. ([5]) is due to the 
extra pinching of the electron beam by the transverse 
electric field E^. We note here that E^ counters the 
magnetic expulsion of the ambient plasma. At the same 
time it reinforces magnetic pinching of the beam. The 
generalized momentum conservation in the z— direction 
gives 



IjVjz = Jjo^jzo H {Az - Azo) . 



FIG. 1: (color online) Snapshots of the the evolution of trans- 

(5) verse electromagnetic Weibel fields ( Ex and Bx) and beam 
filament densities (nt) during the nonlinear stage at a time 
T = 20(27r/a;pe) for four different simulation cases: (a) Cold 
electron beam in a collisionless background plasma, (b) cold 
e-beam in a collisional background plasma and (c) hot elec- 
tron beam in a collisionless background plasma and (d) hot 
electron beam in a collisional background plasma. 

The relativistic treatment of the instability could be 
somewhat cumbersome. However, the essential physics 
can be learnt from the non-relativistic equation of mo- 
tion. For a warm electron beam the equation of motion 

(6) reads as 



If there is no dissipation, then we may derive conservation 
laws for the system. From Eqs. ([2|), ([5]) and (j6|), we have 



dvb± Vbz 
m — 7— = — e — 



dt 



-V, A, 



V±P 

rib 



(8) 



3 



where P is the beam pressure related to the beam emit- 
tance. On hnearizing Eq. ([8]) for smah magnetic field 
perturbation, we have 



'b± 



dt 



rib 



(9) 



where Az = eAz/mc^, Po = {vbz -Vpz)/c^ Vbz/c, and it 
was assumed that \/±5P = 3v^^\/ ±Snb. Eq. ([9|) together 
with the continuity equation dtSub = —rib (V^ • Svb±) 
yields 



(10) 



where = 3v^^ is the square of the beam's sound speed. 
Eqs. ([2]), (|3]), and (p!Q|) form a set of equations to de- 
scribe the sound like perturbation/filamentation of the 
electron beam density. The simplest case of collision- 
less plasma {u = 0) and long wavelength perturbation 
(I k± |<C a;^/c^, cOp being the ambient plasma frequency) 
gives 



^ = 0, 
rib 



(11) 



where = — c^pQUb/riQ is the modified sound speed. 
One may note that the cold beam (c^ < (? jS'^rib / no) is 
unstable due to the Weibel instability whereas the warm 
beam is stable. Under the warm beam approximation 
the dispersion for sound waves is given by cj^ = c^/c^. 
These waves are stable and the Weibel instability does 
not occur for sufficiently high transverse beam temper- 
ature and low beam/plasma density ratios. With finite 
plasma resistivity taken into account the dispersion rela- 
tion for the sound waves reads as 
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where uob is the beam plasma frequency. For large scale 
perturbations {k^ <C /c^g, k~^ = c/uOp) and small colli- 
sion frequency, Eq. ([12]) yields two damped modes and 
one growing mode given by 
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Thus collisions drive negative energy waves in the system, 
leading to the Weibel instability of a warm electron beam. 



which would be stable in collisionless plasmas. This re- 
sult has a profound impact on the understanding of the 
Weibel instability in plasmas. 

To check our analytic theory, we carry out detailed 2D 
PIC simulations. The relativistic electron beam prop- 
agates in the negative z-direction with the initial ve- 
locity Vb^z' The compensating return current of the 
ambient plasma electrons fiows with the initial velocity 
v^. The plasma ions are immobile and have the den- 
sity no = rib ^ rip. The simulation domain has the size 
X xY = (20 As X 20 As), where As = c/uOp^ is the plasma 
skin length. All simulations are performed with 64 par- 
ticles per cell and with a grid size of, 5x = 5y = 0.125 As- 
The density ratio between the beam and plasma electrons 
is Up /rib = 9, whereas the beam and background plasma 
electrons have velocities Vb = 0.9 c and Vp = 0.1c. The 
binary collisions are simulated with a newly implemented 
collision module in the relativistic PIC code Virtual Laser 
Plasma Laboratory (VLPL) [l3]. We record the evolu- 
tion of field energy for every component Fi of the fields 
E and B as {eFirrieCUJpe)'^ dx dy where eFi/rUeCuOp 
represents the relativistic field normalization. 

We use the electron beam initial temperature of T5 ~ 
70 keV and the ambient plasma collision frequency 
^ei/^p = 0-15 for these simulations. 
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FIG. 2: (color online) Time evolution of the perpendicular 
and parallel Weibel E B field energies {E'i, B'i, for 
four different simulation cases as described in the previous 

Fig.m 

Fig.[T]shows snapshots of the transverse E and B ffelds, 
and the structure of the beam ff laments at the time, 
T = 20{27r/ujpe) for four different cases: (a) cold electron 
beam and collisionless background plasma, (b) cold elec- 
tron beam and collisional plasma (c) hot electron beam 
and collisionless plasma and (d) hot electron beam and 
collisional background plasma. The beam density ffla- 
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mentation is shown in the last column in each panel. In 
collisionless case (a), the filaments are small, comparable 
with the background plasma electron skin depth. In the 
collisional case (b), the filament size is bigger. This can 
be explained as a collisional diffusion of plasma electrons 
across the self-generated magnetic fields. In the third 
panel of figure, simulation case (c), the electron beam is 
hot with the transverse temperature T5 = 70 keV, and 
the background plasma electrons are collisionless. Here 
we see no filament formation. The temperature of the 
electron beam stabilizes the Weibel instability. Physi- 
cally the thermal pressure of the electron beam prevails 
over the magnetic pressure in this case. Hence, the mag- 
netic field pinching which actually drives the instability 
does not occur resulting in the suppression of the Weibel 
instability. The last panel of the figure depicts the fil- 
ament formation in the simulation case (d), where elec- 
tron beam is hot and the plasma is collisional. Although 
the beam temperature is the same as in the stable colli- 
sionless case (c), the background plasma collisions revive 
back the instability. This is due to the collision induced 
generation of negative energy waves as discussed earlier 
in the theoretical model. 




FIG. 3: (color online) Beam filaments from 3D PIC simula- 
tions. The ambient plasma is cold and collisionless, while the 
electron beam is hot and collisionless. The other parameters 
are same as in Fig. [B^c). 

Fig. [2] shows the evolution of electric and magnetic field 
energies in the four cases corresponding to the simula- 
tions in Fig. [H The energy axises in Fig. [2] use logarith- 
mic scales. We see a stage of linear instability, where 
the field energies build up exponentially in time. It is 
followed by a nonlinear saturation. The linear instabil- 
ity stage is present in the simulations (a), (b) and (d). 
The simulation (c), where the electron beam had high 
temperature and the background plasma was collision- 
less, shows no linear instability and no significant build 
up of the magnetic field energy. 

The linear growth rates calculated from the simula- 
tions results agree well with the theoretical model. After 
the linear stage of the instability, an electrostatic regime 



of the filamentation instability characterized by the mag- 
netic attraction of the filaments starts and the the field 
energies saturate. Some small fluctuations around the 
saturated field energies can be seen. These fluctuations 
occur due to the collective merging of the filaments as 
also discussed in [9]. 

We also have done a number of 3D PIC simulations of 
Weibel instability, varying the beam temperature and the 
plasma collision frequency. To our surprise, we found no 
stabilization even in the collisionless case for high beam 
temperatures. The corresponding simulation is shown in 
Fig. [3l Although the electron beam in this simulation 
had the high transverse temperature, and the plasma 
had no binary collisions, we see a lot of filamentation 
due to the Weibel instability. We explain this fact in 
terms of anomalous plasma collisionality. Indeed, there 
is an oblique mode in the 3D geometry which couples 
the Weibel and the two-stream instabilities The 
two-stream mode generated electrostatic turbulence in 
the plasma. Stochastic fields associated with this tur- 
bulence scatter the beam and plasma electrons and lead 
to an effective collisionality in the return current. This 
anomalous effect revives back the Weibel instability. It 
may be noted here that we have always taken the back- 
ground plasma as cold. The allowance of finite back- 
ground plasma temperature decreases the growth rate 
of the instability. The interplay of collisions in different 
regimes of beam and background plasma temperatures 
can be found in Refs.[5]. 

In summary, we have developed a simplified model 
which identifies the collisional Weibel instability as the 
instability of the unstable negative energy mode driven 
by the collisions in the background plasma. An impor- 
tant result of this study is that beam temperature does 
not kill the Weibel instability in the presence of collisions 
in beam plasma system. An alternate explanation on the 
persistence of the Weibel instability in 3D geometry is of- 
fered. It is attributed to the anomalous collisionality of 
the beam-plasma system due to the two-stream mode. 
We have also derived the conservation laws for the beam 
plasma system which are useful for benchmarking the 
numerical codes. Detailed 2D simulations on the Weibel 
instability of an electron beam in two-dimensional trans- 
verse geometry have been performed, which essentially 
confirm the theoretical prediction. 
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